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The color suppressed modes B*^ I)^*^''ry^'' are analyzed in perturbative QCD approach. We 
find that the dominant contribution is from the non-factorizable diagrams. The branching ratios 
calculated in our approach for B" D^*'''^rj agree with current experiments. By neglecting the 
gluonic contribution, we predict the branching ratios of _B" D^*^'^r]' are at the comparable size of 
B° D^*^°TT°, but smaller than that of B° D'-'''>°r]. 
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I. INTRODUCTION 



The hadronic decays i?° — *■ D^*'>^ri^''> are color suppressed modes, which belong to class II decays in the factorization 
approach (FA) P| . The relevant effective weak Hamiltonian for these decays is given by 

Kff = % VcbV:Jc^{^,)o^{^,) + c■2{^i)02i^i)\ , (i) 



(N ■ V2 

> , 

where the four-quark operators are 
O ' Oi = {db)v-Aicu)v-A , O2 = {cb)v-Aidu)v-A , (2) 

i> : 

■ with the definition {qiq2)v-A = 9i7m(1 ~ 75)92- The Wilson coefficients Ci ~ —0.2 and C2 ~ 1 are calculated at mi, 
' scale. The main contribution of these decays in FA is proportional to the Wilson coefficients 02 = Ci + C2/3, which 
is a small number. That is the reason why class II decays usually have small branching ratios. Theoretical study 
' of -B° — > D'^*^^r]^'^ decays gives a branching ratio of 10~^ 0- However, recent experiments by Belle and BABAR 
Qh! show that the branching ratios of class II decays are not so small 0, Q. The branching ratios of i?° — > l)^*^°r/ are 
of 10~^. Although the gluonic mechanism can enhance the D^rj' decay branching ratio to 10"'^ 0, it may be 

(]J [ difficult to explain the large branching ratio of J3° D^*^^r]. It means that the non-factorizable contributions in these 
• decays are very important. This is confirmed in the recent theoretical study on charmed final state B meson decays 
^ ' in perturbative QCD approach 
• fH , The perturbative QCD approach (PQCD) for the exclusive hadronic B decays was developed some time ago [3,|1|, 
/\ ' and applied to the semi-leptonic (3 and non-leptonic decays ^3 ^3 successfully. In this formalism, factorizable 
^ , contributions, non-factorizable and annihilation contributions arc all calculable. By including the fcy dependence of 
- - ' the wave functions and Sudakov form factor, this ap pro ach is free of endpoint singularity. Recent study shows that 
PQCD approach works well for charmless B decays [Tol fllllT^ . as well as for channels with one charmed meson in 
the final states Hill. We will show the PQCD calculation of -> D^*^°r)^'^ decays in the next section, and discuss 
the numerical results in section ITlTl The conclusion is in section Hvl 



II. B° ^ D^*^°r)^'^ DECAY AMPLITUDES IN PQCD 



In two-body hadronic B decays, the two outgoing mesons are energetic. Each of the valence quarks inside these 
mesons carries large momentum. Most of the energy comes from the heavy b quark decay in quark level. The light 
quark (d quark) inside B^ meson, which is usually called spectator quark, carries small momentum at order of Aqcd- 
This quark also goes into final state meson in spectator diagrams. Therefore, we need an energetic gluon to connect 
this quark to the four quark operator involved in the b quark decay. Such that the spectator quark get energy from the 
four quark operator to form a fast moving light meson. The hard four quark dynamic together with the spectator quark 
becomes six-quark effective interaction. Since six-quark interaction is hard dynamics, it is perturbatively calculable. 
The non-pcrturbative dynamics in this process is described by the wave functions of mesons consisting of quark and 
anti-quark pair. The decay amplitude is then expressed as 

Amplitude - [d^kid''k2d^k3 Tr[C(t)$s(fci)$c(.) (fc2)$^(') (^3)1^(^1, ^2, fca, i)e-^(*)] . (3) 
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FIG. 1: Color-suppressed emission diagrams contributing to the _D'-*'''r;'-'' decays. 



Here C{t) is the QCD corrected Wilson coefficient of the relevant four quark operator at scale t. Although next-to 
leading order results have been given we will use leading order one here ^J. $i are the meson wave functions, 
which include the non-perturbative contributions in these decays. The non-perturbative wave functions are not 
calculable in principal. But they are universal for all the hadronic decays. We will use the ones determined from other 
measured decay channels |a lift IllL . The exponential S{t) is the so-called Sudakov form factor, which includes the 
double logarithm resulting from the resummation of soft and collinear divergence. This form factor is also calculated 
to next-to lea ding order in the literature The Sudakov factor effectively suppresses the soft contributions in the 
process [allOllllj. thus it makes the perturbative calculation of hard part reliable. 

Now the only left part of the decay amplitude is the hard part H{t). Since it involves the four quark operator and 
the spectator quark connected by a hard gluon, it is channel dependent, but perturbatively calculable. There are 
altogether 8 kinds of diagrams in our BO £)(*)077(') decays, which are shown in Fig. 1 for the spectator diagrams and 
Fig. 2 for the annihilation type diagrams. Notice that in Fig.l, the yy*^'-* meson consists of dd content, while in Fig. 2, it 
is a uu pair making r]^''> meson. Since ry^'^ meson is isospin singlet {uu + dd), these two sets of diagrams give relatively 
positive contributions. On the other hand, in case of — D^*'>°tt° decays 0, where 7r° is isospin triplet (uu — dd), 
these two sets of diagrams give destructive contributions there. Fortunately, as we will see later in the next section, 
the annihilation type diagrams are suppressed comparing to the spectator diagrams. Therefore, the branching ratios 
of these two kinds of decays are still comparable. 

The structures of the meson wave functions are 



v'oLiP) 



B,n{P) ■■ + mBh5<pBix) , 

DoutiP) ■■ isi^ + mD](l)D{x) , 
Dl^tiP) ■./[^ + mDA^D'{x) , 
Isi/P't'Aix) + ■mo(f>p{x) + Ctoo(A- A+ - ^)4'Tix)] 



(4) 
(5) 
(6) 

(7) 



with Too = TOj/(TOtj -|- nid) = 1.4 GeV, utilizing isospin symmetry. And the light-like vectors are defined as n+ = 
(1,0, Ot) and n_ = (0, 1,0t). As shown in ref.Q, (t^B is identified as 0+ and the contribution of another B meson 
wave function (j)B oc ~ 'P- is smaller in the PQCD calculations, therefore we neglect it. Applying for heavy quark 
symmetry, there is only one independent distribution amplitude (jy^i*) in the heavy D^*^ meson wave function [a Il3l|. 
However, there are three distribution amplitudes for the light -q^'^ meson wave functions like the tt meson wave 



-1 arc for ?7^2t 



with u {d) carrying the momentum X3P3, while C = — 1 for 77^2* with 



function. The coefficients ( = 
(d) carrying the momentum X3P3. 

The gluonic mechanism of 7/ may make sizable contributions in the B D^*^!^' decays, but they are usually model 
dependent Thus we will not consider it here. The and 77' mesons are mixtures of fiavor SU(3) octet (jyg) and 
singlet (770) states in a two-mixing- angle formalism ^6^, 



rj = cos 6^81778) - smfc'ol'yo) , 
77' = sin 6*817/8) + cos 6*017/0) ■ 
The definitions of the decay constants of 77 and 77' are as follows: 

{Q\u-t^-t,u\ri^'\p)) = if;,,p^, {Q\d-t^-t,d\r^^'\p)) = z^.p^. 

The ss components of rj and 77' are not relevant in our decay channels. Therefore we did not show them. The decay 
constants in the two-angle mixing formalism are 



(8) 
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/" = /,f = Acos08"^sin0o, 
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FIG. 2: Annihilation Diagrams contributing to the — + D^*^'^rj^'^ decays. 




The parameters are determined to be Qi 



98 = -22° ~ -21°, fs = 1.28A, 

9o = _9°~-4°, /o = (1.20 - 1.25)A . 



The 7? and rj' meson can also be expressed as a mixing in the quark flavor basis [l7| 



cos a — sm a 
sin a cos a 



(wM + dd)/^/2 \ 
ss J ' 



(12) 



(13) 



where a ^ n + 6 — arctan(l/V2) describe the deviation from ideal mixing. The angle 9 is one-angle mixing parameter. 
From the above eqn. (|13|l . applying isospin symmetry, we have 



/,; = /,f-/xC0sa/V2, 



The range of mixing parameters are determined to be 9 
The B° — > D^*^^r]'--'^ decay rate has the expression as, 



-17° — 11° Ei 
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Including the hard part and the meson wave functions, the — > Z?^*)^?]^') decay amplitude is written as 

M{B° ^ 5(*)°r/')) - /c(,)6nt + /bCcxc + Mint + Xcxc , 



(14) 
(15) 



(16) 
(17) 



where Jb = 190 MeV, fo ~ fo' = 240 MeV are the B and D^*^ meson decay constants, respectively. The functions 
^i„t, and ^cxc denote the internal M^-emission, and VF-exchange contributions, which come from Figs. 1(a) and 1(b), 
Figs. 2(a) and 2(b), respectively. The functions A^int, and Aioxc represent the internal H^-emission, and VF-exchange 
contributions, which come from Figs. 1(c) and 1(d), Figs. 2(c) and 2(d), respectively. The expressions of the four 
functions are already shown in the appendix of ref. ^ for B Dn decays. One need only replace the pion wave 
function by the rj^''^ wave function in those expressions. 

In FA, only factorizable contribution of ^i„t (Fig.l(a)(b)) has been considered. Since ^i„t is proportional to the 
small Wilson coefficient 02 = Ci + C2/3, the branching ratios predicted in FA is smaller than the experiments. Now 
in PQCD approach, all the topologies, including both factorizable and nonfactorizable ones, and also annihilation 
type ones have been taken into account. In fact the non- factorizable contribution Mint, which is proportional to the 
large Wilson coefficient C2/3 is the dominant contribution in the B^ — > D'^*^^!]'''^ decays. The reason is that the two 
non-factorizable diagrams in Fig. 1(c) and (d) do not cancel each other like the B to two light meson decays, where 
the distribution amplitudes of wave function are symmetric [lol lll| . The large difference of c and u quark mass makes 
the contribution of Mint large. Very recently, the soft coUinear effective theory also confirms that the non-factorizable 
Mint dominate over the contribution of f]j(t)£,int |lS| . 



III. NUMERICAL RESULTS 



As stated in the above section, we need various wave functions in our numerical calculations. Considering the 
previously calculations of other decay channels 0, 0, ITill Ull IT^ , the B meson wave function has been determined as 



4>b{x, b) = Nbx'^{1 — x)"^ exp 



xMe 



OJB 



,2 



(18) 
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TABLE I: PQCD predictions with one angle mixing formalism (I) and two angle mixing formalism (II) and experimental data 
(in units of 10^") of the B° branching ratios. 



Decay mode 


PQCD (I) 


PQCD (II) 


Belle 


BABAR 


PDG 




1.7 ~ 2.3 


2.2 - 


- 2.6 






< 9.4 


B" ^ &'ri 


2.4 ~ 3.0 


2.6 - 


- 3.2 




2.41 ±0.50 




B" ^ D^n" 


2.3 ±0.1 






3.1 ±0.6 


2.89 ±0.48 






2.0 ~ 2.7 


2.6 - 


- 3.2 






< 14 


B° ^ D'"ri 


2.8 ~ 3.5 


3.1 


- 3.8 










2.8 ±0.1 






2.7 ±0.9 







where the shape parameter is chosen as = 0.4 GeV. The normalization constant A''^ is related to the decay constant 
Jb through 

j dx4>B{x,Q)^^. (19) 

The meson distribution amplitude is given by 

(x) = -^/b(.,x(1 - a;)[l ± (1 - 2a;)] , (20) 

with the shape parameter Cd — Cd* — 0.8 ± 0.2 The range of C^ic.) was extracted from the B D^*Hv decay 
spectrum at large recoil assuming lob — 0.4 GeV for the B meson wave function 9]. We do not consider the variation 
of with the impact parameter 6, since the current data are not yet sufficient to control this dependence. The 

light ry*-'' meson wave functions are chosen to be the same as the pion wave function according to isospin symmetry, 
since the relevant valence quarks here are mainly uu and dd. 
In the numerical analysis we adopt 

A^"*^ = 250 MeV, Mb = 5.2792 GeV, Mw = 80.41 GeV. 

Choosing \Vcb\ = 0.043 and |Kd| = 0.974, we obtain the PQCD predictions for the B° branching 
ratios shown in Table For comparison, we also list the 5° D^*^^Tr° decay branching ratios in this table. The 
theoretical uncertainty comes from the variation of the shape parameter for the meson distribution amplitude, 

0.6 < Cb>(,) < 1.0 and the t] and rj' mixing parameter 9 — —17° ~ —11° for the one-angle mixing formalism (I). 
The range of parameters of two- angle formalism (II) are shown in ea. (|12|l . From numerical study, we notice that the 
branching ratios do not vary much upon the variation of Cjj(,) . This can be seen from the numbers of B^ — > Z)(*)°7r*' 
m Tabled since it only depends on this parameter. Most of the uncertainty of — > D'^*^'^ •q'^'^ decays is from the 
mixing parameter 9 or 9^ and ^o- The branching ratios of — > D'^*^^r]' increase while B^ — > l)'*)"/? decrease as the 
angle 9 getting larger. Other input parameters, such as parameters of B meson wave function also affect the branching 
ratios, but they are mostly constrained by other well measured decay channels, like B ^ rnr |ll,] and B Kn 10] 
decays, etc. The uncertainties of PQCD approach itself mainly come from the unknown higher twist contributions 
and higher order calculations of corrections. No numerical estimation of higher twist contribution exist, although 
it is expected to be suppressed. The higher order QCD calculation in i? — > (pK decay shows that next-to leading order 
as correction may not be small in certain channels 19]. 

The recently observed class-II decay B^ D'^*^^ri branching ratios are also listed in Table |T| H S H^. It is easy 
to see that, our results agree with the experimental measurements within errors. Since we do not consider the extra 
gluon fusion contribution to B^ D'^*^^rj' decay, the not yet measured B^ D^*^^r]' branching ratios are a little 
smaller than the B'^ D^*^°r] branching ratios. But they are still comparable with the B^ l)(*)°7r*' branching 
ratios. The reason for this comparable result is that we apply the assumption of exact isospin symmetry. We use the 
same wave function for yy^'-* and tt meson where the only difference is the decay constant. The difference of dynamics is 
the constructive or destructive contribution from annihilation type diagrams. This destructive contribution makes the 
B^ l)(*)°7r'' branching ratios smaller than that of B^ D'-*^^r] decays. The numerical results also show that the 
dominant contribution comes from the non-factorizable contribution Mint- The factorizable contribution foS^int and 
annihilation contribution M^xc are only 20-30% of Mint- Factorizable annihilation contribution fs^exc is negligible. 
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IV. SUMMARY 

In this work, we calculate the branching ratios of D'^*^'^ rj'^'^ decays in the perturbative QCD approach with fcy 

factorization, which is free of endpoint singularity. Belonging to class II decays in FA, these decays receive dominant 
contributions from the non-factorizable diagrams. Naive factorization breaks down in these color suppressed modes. 
The branching ratios calculated in our approach for — > D^*^^r] agree with current experiments. We predict the 
branching ratios of B^ D'^*^^ri' without gluonic contributions are at the comparable size of B'^ D'^*^^tt^ , but 
smaller than B" D'-*^°r]. They may be measured soon in the B factories. 
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